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A method of construction of error-correcting codes is described as well as a class of linear 
q-ary codes. 

1, Introduction 

This paper generalizes the results set forth by the author in [1]. At the beginning we describe a meth- 
od of construction of error-correcting codes that makes it possible to interpret a linear q-ary code of weight 
d as a set of rational functions of degree >d. The possibilities of the method are illustrated by the construe- 
tion of a wide class of codes called (L,g)-codes. The codes have the following parameters: n < q m , k >n 
— 2mt, and d > 2t + 1. The decoding of these codes is not more complicated than the decoding of Bose — Chaud- 
huri-Hocquenghem (BCH) codes; for some (L,g)-codes there exist also simpler decoding schemes. In con- 
clusion we shall consider some particular cases of (L,g)-codes, and we shall study their behaviour for n 

— ► OO # 


2. Rational Representation of Codes 

1. Principal Isomorphism. Let q = p^, p is a prime, L = {a^, a 2 , . . . , oi n } , n < q m , ct[ 6 GF(q m ), 

^ Oj, with GF(q in ), being the minimal field containing L, S is a vector space over GF(q) of dimension n, Sjj 
is a metric space S with a Hamming norm, and R a vector space of rational functions over GF(q m ) of the 
form 

E(*)= i ^GF(q),a^L. (1). 

i=l * 

Let us introduce a norm in R as follows: for any irreducible fraction 

= t + Q, |l(z)|| = degq>( Z ). (2) 

\Z) 

By R^ we shall denote a space R that has been metricized with the aid of this norm . Let x 6 S, x = (bj , b 2 ...» 

b n>- 

The mapping x-+l, x (z) = 2^i/( z — a i) is an isomorphism of S onto R. Since £ x (z) =11x11 , this mapping 

i=l 

will be an isometry of Sjj onto R d , i.e., S H « R d . This isomorphism makes it possible to define codes as 
subsets of R. In particular, a linear code of weight d is a linear subspace of R consisting of fractions of 
degree not smaller than d. 

2. Characterization of R for a Prime q. If q is a prime, then GF(q) will coincide with the residue 
field Z~of integers modulo-q, and it will be useful to characterize the elements of R as follows. 

To each vector x 6 S, x = (b*, b 2 , . . . , b n ) we assign a polynomial f x (z) = (z-a^b^z-a^k * 2 3 . . . (z— o; n )kn. 
Then any fraction £ 6R can be represented in the form £ = f x '/f x , where x 6 S and f x ’ is the formal derivative 
of the polynomial f x . 

3. Construction of Binary Codes. For the construction of binary codes it is possible to use a more 
general procedure. Let S be a vector space of dimension n over GF(2), K an integer domain with unity, 
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L a set of n distinct elements K: L = {a v a 2 , . . . , a n ] , D a multiplicative monoid of polynomials over K with 
roots in L, and D 1 2 a submonoid of D consisting of polynomials with roots of even multiplicity. 

Two polynomials in D that differ by a factor belonging to D 2 will be assumed equivalent. A factor 
monoid D/D 2 with respect to this equivalence relation is a multiplicative group. As representatives of 
equivalence classes we shall take polynomials of smallest degree in a class. Such a representative can be 
obtained from any polynomial f 6D by reduction modulo-2 of all the multiplicities of the roots of this poly- 
nomial. 

Let x 6 S. The correspondnece x — * f x (z) is an isomorphism of the additive group S onto the multipli- 
cative group D/D 2 : S^D/D 2 . Here degf x (z) = d if the weight of x is equal to d. By using this isomorphism, 
it is possible to define a linear binary code of weight d as a multiplicative subgroup of the group D/D 2 con- 
sisting of polynomials of degree not smaller than d. 

For obtaining such a code, it suffices to specify a monoid F such that D 2 cz F cz: D. Then the linear 
code C will be equal to F/D 2 . Let us point out the wide possibilities of such a method of code construction. 
Binary codes are defined with the aid of polynomials over an arbitrary integer domain (in particular, it is 
possible to use real polynomials) , and the weight of the code vector is determined on the basis of the degree 
of the polynomial. 

By taking K = GF(2 m ), we obtain f'(z) = 0 for any f € D 2 . Since for any <p, ip 6 D/D 2 we have (<; pip )'/<p$ 

= <p'/<p +$'/$, it follows that in this case D/D 2 will be isomorphic to an additive group R of functions of the 
form <p'/<P> and we obtain the principal isomorphism. 

3. (L , g) -Codes 

1. Definition. Estimation of Weight and Power. We shall use the principal isomorphism for defining 
a wide class of linear codes. Let g(z) be a polynomial over GF(q m ) that has no roots in L. For any | (z) 

= ip(z)/cp(z) £R the polynomial <p(z) is relatively prime to g(z), and hence invertible in the algebra G of poly- 
nomials modulo-g. Let us define an (L,g)-code as a set of elements ^ (z) 6R such that £ (z) — Omodg(z). 

The linearity of such a code is evident. Estimates of the weight and power are presented in the fol- 
lowing theorems. 

THEOREM 1. An (L,g)-code has not more than m -deg g(z) check symbols. 

Proof. By definition, an (L,g)-code is the kernel of a linear mapping of the space R into an algebra 
G of polynomials modulo g(z). Thus the number of check symbols coincides with the dimension over GF(q) 
of the image of R under this mapping, and this dimension cannot be larger than m - deg g (the dimension of 
the algebra). 

THEOREM 2. The weight of any element of an (L,g)-code is not smaller than deg g + 1. 

Proof . If i = 0 mod g and £ = ip/f, then ip = Omodg, i.e., degip >degg, and hence deg£ >degg + 1, which 
completes the proof. 

For q ^ 2, (L,g)-codes will be more effective that have a check matrix with an additional row of ones 
(modified (L, g)-codes) . 

THEOREM 3. The parameters of a modified (L,g)-code satisfy the conditions n < q m , k > n-(2t-l)m 
-1, d > 2t + 1. 

n 

Proof. Let x6S, x = (b, , b 2 , . . . , b n ) be a code vector. Then £ = i p/cp = Omogg, S b i = °- Since the 

leading coefficient of the polynomial ip is 2 b i> lt follows that deg cp > tp + 2 >deg g + 2. 

1 

Hence for the correction of t errors it suffices to take a polynomial g such that deg g = 2t-l. 

2. The Check Matrix. In connection with (L,g)-codes there appears a new form of check matrix, i.e., 
the check matrix is in the form of a row of polynomials (a check matrix of compact form). Such a form is 
convenient in the analysis of the properties of a code and in the construction of decoding schemes. The check 
matrix of (L,g)-codes can be represented in the following forms. 
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First compact form: 


7\c = 


second compact form: 

Here r = deg g(z). 
Expanded form: 


1( } g(z)-g(a n ) 
z — cq v z — ct n v n 


^2C — 


z r - ar z r 

L o 1 1 


-g 




f " 1 («n) ■ 


T = 


g- 1 K) • ■ • • 


g -1 («i>ai- • ■ 

• -g^( a«)a» 

g~ l K) aj- 1 • - 

• ■ S" 1 («») a n 1 


We can convince ourselves of the validity of these forms in the same way as in [1], 

Thus the check matrix of an (L,g)-code can be obtained from a "rectangular Vandermonde matrix" by 
multiplying the columns by the elements g' 1 ^). It is possible to define more general codes by taking as 
factors any nonzero elements not necessarily equal to g -1 (^i). As before, any r columns of the thus-obtained 
matrix will be linearly independent over GF(q in ); therefore the weight of the code will be larger than r. 

By taking as factors the elements g -1 (cq), we obtain the definition of (L,g)-codes given at the beginning 
of § 3, and we can estimate the parameters with the aid of more refined methods, not related to the check 
matrix (Theorems 1 and 2). Thus if q = 2 and g does not have multiple roots, then the weight of any code 
element will be larger than 2r; in this case the check matrix Tp will yield a bound that is larger than r only. 

3. Decoding. Let y = x + e, with y, x, e 6 S, and let x be the transmitted code vector and e the error 
vector with nonzero coordinates bi 1? bi 2 , . . . , bt^, k < t. By multiplying the vector y by the matrix T^ (the 

k k 

transpose of T^), we obtain the syndrome S(z) = 2 b iZ [(g(z)-g(^iz)]/ (z-Qiz) Let ( z 

i=i 

-aij). Then 

fS = x \> mod g, deg / ^ Z, deg o|) < Z, ( 3 ) 

and for decoding it suffices to solve the congruence (3) for f and ip. 

THEOREM 4. If deg g = 2t, then for any polynomoal S(z) of degree < 2t the set D of pairs (u,v) such 
that uS= vmod g, deg u < t, deg v < t will contain a nonzero element. All the pairs (u,v) belonging to D 
will have the same ratio v/u = 77 . 

proof. Let A be a linear operator of multiplication by S in the coset algebra of polynomials modulo 
g(z). In the basis 1, z, . . . , z 2t-1 the matrix of this operator has the form 

2t — i 

, where z y S ( z ) = ^ c n z% * 7 — 0* L •••> 2t — 1. 
i=0 

Let us denote by At ( t+i the submatrix of A constructed from t upper rows and the first t + 1 columns. The 

equation ’ ■ w = o < 4 > 

always has a solution, whence follows the existence of a quanitity u that satisfies the condition of the theorem 
and also its method of determination. 

Now let (UpV^, (u 2 ,v 2 ) €D. Then UjS = vqmodg, v 2 = Su 2 modg, where the polynomials S, v l5 v 2 and g 
were obtained from S, v 1# v 2 and g by division by the largest common divisor of S and g. Since S and g are 
relatively prime, it follows that u t v 2 = v 1 u 2 modg, v t u 2 = Omodg, and since the degree of the poly- 
nomial in the left-hand side is smaller than the degree of g, it follows that VjAq = v 2 /u 2 , so that Vj/uj = v 2 /u 2 , 
which completes the proof. 

From this theorem we obtain the following rule of decoding of (L,g)-codes; 1) Find a syndrome S(z); 
2) find a solution u of system (4); 3) calculate v = u ■ Smodg; 4) find 77 = v/u (this ratio will coincide in 
the case under consideration with ip/f from the congruence (3), since f and ip are relatively prime); 5) de- 
compose 77 into simple fractions (for this purpose it is necessary to find the roots of f(z) and calculate c\i 
= ip (ail) /i' (aii)); 6) determine the values of the errors: biz ~ ~ c iZ- 


C 2f-1 0 • ■ 

* • C 2t-1 Zt-1 

coo ; . 

• C 0 2/-1 
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Another method of decoding involves the use of the check matrix T 2C . By multiplying the vector y by 
the matrix T 2C T we obtain the syndrome S (z) = 2 b iZ (z r — aj^/tz-a^) Let 

■f— 

Then the decoding reduces to solving the congruence 

f-S = ^ modz r (®) 


for ail ^ 0,1 = 1>2» ...» k, or the congruence 
for = 0 for some j. 


f-S = i|)mod2 r_1 


( 6 ) 


In the first case deg f = k, and in the second case deg f = k— 1, so that for r = 2t there exists by virtue 
of Theorem 4 a unique solution (f, ip) of the congruence (5) with f and ip being relatively prime, deg f < t, deg 
ip < t, or of the congruence (6) with f and ip being relatively prime, deg f < t — 1 , deg ip < t— 1. 


Berlekamp (2) has proposed an iterative decoding procedure for BCH codes that is in fact a method of 
solution of the congruence (5). The procedure involves successive transition from a congruence modulo z 2 
to a congruence modulo z 3 , then to a congruence modulo z 4 , etc. The Berlekamp algorithm saves a large 
amount of equipment, since it does not require storage of the matrix Aj. f. +1 . 


5. Special (L,g)-Codes. 

lows: 


Below we shall consider some particular (L,g)-codes that are related as fol- 


i 

Cumulative 

codes 

i 

BCH 

codes 


-{L, 5)- codes ^ 

Separable 

codes 

1 

i i 

Srivastava codes | Irreducible 

i | codes 

Gabidulin codes 1 


The classification is based on the form of the generating polynomial. Cumulative codes are of the 
same type as BCH codes. Separable codes have a peculiar decoding scheme. For both these classes of 
codes it is possible to improve the parameter bounds in the binary case: n < 2 m , k > n— mt, d > 2t + 1. 


The Srivastava codes and the Gabidulin codes were obtained in 1967. Irreducible codes have the 
"strongest" algebraic structure, and we have the impression that in studying the properties of these codes 
it is possible to make good progress. 


4. Cumulative Codes 

Cumulative codes are codes whose generating polynomial has one root: g(z) = (z— a) r . 

The maximal!/ that can be taken for such codes corresponds to GF(q m ).— {a} , so that their maximal 
length is equal to qi* 1 — 1. A particular case of such codes for g(z) = z r is BCH codes (with such a choice 
of g(z) the parity check matrix T e goes over into the well-known parity check matrix for a BCH code). 


THEOREM 5. All the cumulative codes with the same r have the same spectrum. 

Proof. If x is an element of a (GF(q in )-{a} , (z-a;) r )-code, then £ x (z) = ip(z)/<p(z) = Omod(z-a) r , so 

that 

(*+<*) /V , r 

-^'Vaj'=° m0dZ ' 

If <p(z) has roots aij, . . . , aij, then cp ( z + a) will have roots aij— q, . . . , aii—a that lie in GF(qm) — {o}. 
Thus there exists aweight-preserving one-to-one correspondence between any cumulative code and a BCH 
code. 


THEOREM 6. BCH codes are the only cyclic (L,g)-codes with L = {\,a,o? , . . . , a 11-1 }, where a is the 
n-th primitive root of unity. 

Proof. With such a choice of L it is possible to characterize cyclic codes with the aid of a rational 
representation as follows: a cyclic code is a set of fractions in R (§ 2) that form a linear space, this set 
being closed under the substitution z -► az, a GL: C = {£ (z) €R: £(o;z) 6C, a 61/}. We shall assume that the 
generating polynomial g(z) of a cyclic (L,g)-code has a nonzero root p. Any £ (z) belonging to this code sat- 
isfies the congruence 

^ olz ) := '|(S^^ 0modg (2) 
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for all a 6 L; this signifies that ip( z) has in addition to the root p all the other roots of the form a/3, i.e., 
deg^(z) >n, which is impossible. Hence g(z) = z r . 


5. S eparable Codes 

1. Definition. Check Matrix. Separable codes are codes whose generating polynomial does not have 
multiple roots: 

g(z) = (2— Zi)(Z — Z 2 ) ... (Z — Z r ). 


Binary separable codes, just as BCH codes, admit an improvement in their parameter bounds. Indeed, 
in accordance with § 2.2 we have f'/f = Omodg, and since the derivative in a field of characteristic 2 is a 
square, it follows that f f /f = Omod g 2 . Hence for deg g = t we have n < 2 m , k >n — mt, d > 2t + 1. Separable 
codes have the following specific form of the check matrix (Cauchy matrix): T c = II (z^— aj)" 1 II , i = 1, 2, ...» 
r, j = l,2, ...,n, where Z[ and aj are distinct elements of the field GF(q m ) or of an extension of this field. 
By using the check matric T c we obtain an efficient decoding method for such codes based on rational inter- 
polation. 


2. Decoding. Let y, x and e be the same as in § 3.3. By multiplying the vector y by the matrix T c ',we 
obtain a syndrome (£ x , £ 2 > * * • > £ r) , where £ [ = £ (z [) , 


6(2) = 


k 


s 


b,. 


(7) 


In this case decoding involves the reconstruction of the fraction £ (z) on the basis of its values at certain 
points (rational interpolation [3]). 


Let K be a field and £ (z) = i/>(z)/f(z) an irreducible fraction over the field K, deg f = n, deg ip ^ m. The 
order of a fraction is defined as a quantity r(£) equal to 2n for m < n, and to 2n— 1 for m > n. Any irreduc- 
ible rational fraction of order < k can be uniquely reconstructed on the basis of its values at the (k + l)-th 
point of the field K. 

THEOREM 7. Let z t , z 2 , . . . , z k+1 be distinct elements of the field K. Any irreducible fraction ip( z) 
/f(z) = £ (z) 6 K(z) of order r < k can be represented in the form 


£ (z) + 



(8) 




Proof. Among the elements z 1( z 2 , . . . , z k+1 there exists an element z such that f(z[/) 5* 0. Let us 
write Xj = £ (zij). Then 

/ (*) 


*«=*! + 


£1 ( 2 ) 


where £1 (2) = 


^1 (z) ’ ^ z 


If m < n, then r(£) = 2n, degtp < n, r(£ t ) = 2n— 1. If m > n, then r(£) = 2m- 1 , deg^ = m — 1, r(£j) = 2m — 2. 
Hence we have in any case r(£ t ) = r(£) — 1. By applying the previous reasoning to the fraction £ t (z), then to 
the fraction £ 2 (z), etc., we obtain the decomposition (8) after precisely r steps. 

The quantities X 1 ,X 2 » . . . . X r in (8) are determined by successive substitution z = z^, z = z[ 2 , etc. It is 
convenient to perform the calculations by the scheme 


Here 


*1 

y? 




2 2 

yf 

y ( f 



Z Z 

yf> 


y[ 3) 


H 


»?> 

yf 

y?> 




yf> 

y? 



Vi — 

i ( z i-i 

), 



z i+j-l 

. - Zi- 

-1 


yy' — 


-v? 

-l) ’ 


(9) 


2/J 6) 
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If = yii) for some i ^ 0, then the (i + l)-th row must be eliminated from (9). The quantities \ t , 
X 2 , . . . , Ar form the upper diagonal in (9), i.e., Aj = Instead of (9) it is possible to construct inverse 

differences and then use them for determining X y. 


Here p 1 (z 2 z i) = z 2 ~ z i/V 2 “Yi» 


^ y i 
z 2 y 2 
2 S ys 
2 4 y 4 


pi(z 2 Zi) 
Pi(z 3 2 2 ) 
Pi ( 24 ^ 3 ) 


P2 (z 3 z 2 zi) 
P2 (Z4Z322) 


p 3 (Z4Z3Z2Z1) 


P2(2 3 2 2 2l) = 


Z 3 Z 1 

Pi ( 2 3 z a) Pi ( z 2 z i) 


2/2 etc. 


( 10 ) 


In terms of inverse differences the quantities Aj are defined as follows: 

Aj = p j (Zj+iZj . . . 2i) — Pi -2 (2j_ iZj —2 . . . Zi) 

for j > 2 and \ x = y t , X 2 = pj (z 2 z 1 ). The inverse differences Pj( z j+i ■ • ■ z j) are symmetrical in all the variables. 

If q ^ 2, deg g = 2t, k < t, then the fraction £ (z) (7) will have an order < 2t and it can be uniquely re- 
constructed on the basis of its values at 2t + 1 points. If q = 2, deg g - t and k < t, we must effect the fol- 
lowing transformations: 

f(z)=U* *(z)+ZV*{z), V 2 =f, /(* t )gi = f( 2 i ), 

U 2 {Zi)li = V 2 (Zi) (1 +Zili). 

If = 0, then Zi will be a root of v 2 (z), and if 1 + z^f = 0, then ^ will be a root of u 2 (z). Knowing the values 
of yi = u(zi)/v(z i ) = Vzi + | -1 at points z[ at which * 0, it is possible to uniquely reconstruct u(z)/v(z) by 
representing it in the form of a continued fraction (8). 

3. Srivastava Codes. These are codes whose generating polynomial decomposes in a minimal field 
containing L. Such codes are described in [2], where their decoding scheme (based on rational interpolation) 
is also mentioned. Srivastava codes have the following bounds: 

n ^ q m — 2t, k^n — 2 mt, d ^ 2t + 1. 

4. Gabidulin Codes. A particular case of Srivastava codes (for m = 1) was independently obtained by 
Gabidulin [4] who also noted that if Cauchy's matrix (in the case that L c= GF (g) ) is augmented by the unit 
matrix, the code weight will not decrease. Maximal codes defined by such a parity check matrix yield a 
higher rate of transmission than Reed — Solomon codes, i.e. n = q, k = n— 2t, d = 2t + 1. 

5. Irreducible Codes. A separable code is said to be irreducible if the polynomial g(z) is irreducible 

over a minimal field containing L . 

The parity check matrix of an irreducible code consists of the row 


Tu = ( (Zo — Oi) (z 0 — a 2 ) - 1 . . . (Zq — a„) _1 ) , 


where Z 0 is a root of g(z). The value of a rational function over GF(q m ) at a point z 0 6GF(q«ir) determines 
the value of this function at points conjugate to z 0 over the field GF(q m ): 

i(z 0 ) = y Q , g(oz 0 ) = W>,.. • • * l(ff r_l 2 o) =o r -*y 0 . 

All the inverse differences of the same order, calculated for the points z 0 ,c7z 0 , . . . , a r " 1 z 0 , will be conjugate: 


*0 Pi 

ozq <yy o Q-,p i 

a 2 z 0 cr 2 y 0 a 2 pj ap 2 

o 3 z 0 a 3 j/ 0 


P3 


Hence for the decoding of irreducible codes it is not necessary to construct the entire table of inverse dif- 
ferences. It is apparently simpler to perform stage-by-stage decoding of such codes. The number of er- 
rors for a syndrome y 0 can be determined by the following algorithm (here a is the generatrix of the Galois 
group of the field GF(q mr ) over the field GF(q m ) and the notation A: = B has the same meaning as in ALGOL, 
i.e., assign to A the value B): For q ^ 2: 1° write A: = z 0 ; C: = 0; D: = y 0 ; N: = 1; 2° calculate E: = aD 

-D; 3° if E: = 0, go to 5°; 4° write B: = crB; E: = (B-A)/(E + C); C: = aD; D: = E; N: = N + 1; go to 2°; 5° 

print the value of the number of errors equal to N/2; for q = 2: 1° write A: = B: = z 0 ; C = 0; D: = y 0 ; N: 

= l; 2° calculate D: Vd~ ! + A; 3° calculate E: = crD-D; 4° if E - 0, go to 6°; 5° write B: = aB; E: = (B-A) 

/(E + C); C: = crD; D: = E; N: = N + 1; go to 3°; 6° the number of errors is equal to N. 
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In realizing this algorithm it is possible to perform the calculations either in the field GF(q mr )or in 
the field GF(q m ). For running the algorithm it is necessary to perform multiplication in a finite field, then 
obtain the inverse quantity and perform the operation ax = xQ 1 * 1 . 

6. Behavior of (L, g)-Codes for n <*> 

In conclusion let us study the behavior of (L,g)-codes when n — 00 and the transmission rate k/n re- 
mains unchanged. Codes such that d/n — c * 0 are called in this case "good’’ codes, whereas codes with 
d/n— 0 are "bad." It is well known [5,6] that BCH codes, as well as all the codes that are invariant under 
an affine group of transformations, are "bad," whereas "good" codes can be obtained (as was shown by 
Zyablov [7]) by using the concatenation principle discovered by Forney [8]. (L,g)-codes also belongto "good" 
codes. For n — 00 , most of them lie arbitrarily close to the Varshamov-Gilbert bound. 

THEOREM 8. Let L = GF(q m ), n = q m , and let H(x) be the entropy. For any 0 < \ < 1 and e > 0 the 
probability that a randomly selected polynomial g(z) 6L (z) of degree [\n/logn] that has no roots in L will 
generate an (L,g)-code with parameters k/n > 1 — H(d/n) — e will tend to 1 with increasing n. 

Proof . Let us denote by M r a sphere of radius r in the space R^ (§ 2), i.e., a set of fractions in R of 
degree < r. Let Nr be the set of all the numerators of fractions in M r , let V r be the set of all normalized 
polynomials in L(z) of degree t = [An/logn] that do not have roots in L and are divisors of any polynomial 
in N r , and let K be the set of all normalized polynomials in L(z) of degree t that do not have roots in L. By 
A we shall denote as usual the power of the set A. 

Any polynomial belonging to K generates an (L,g)- code C = {£ (z) € R : i (z) = 0 modg (z)} with a trans- 
mission rate k/n > 1 — mt/n > 1 — A., and all g(z) GV r generate codes of weight < r, whereas all g(z) £K— V r 
generate codes of weight > r. Let us estimate the power of all these sets: 

r 

M r = 2 C* B (? - l) 1 ; N r < M r ; V t< N t . C* r ; 
i= 0 

t 

T = 2 c n(— 

i=0 

(as usual, X(n) ~ Y(n) signifies that X(n)/Y(n) — 1). 

Let us determine ju by the condition H(ju) = A.— e, 0 < ju < 1/2 and write r = jxn. All the g(z) €K— V^ n 
generate (L,g)-codes with parameters H(d/n) >H(ju) = A. — 1— k/n— &, and since 

< C £n l0S ”• ^ 0, 

this completes the proof of the theorem. 

Let us note that this result holds also for a narrower class of (L,g)-codes, namely for irreducible 
codes. 
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